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The standard model for molecular recognition of an odorant is that receptor sites discriminate by molecular 
geometry as evidenced that two chiral molecules may smell very differently. However, recent studies of 
isotopically labeled olfactants indicate that there may be a molecular vibration-sensing component to olfactory 
reception, specifically in the spectral region around 2300 cm~^. Here we present a donor- bridge- accept or 
model for olfaction which attempts to explain this effect. Our model, based upon accurate quantum chemical 
calculations of the olfactant (bridge) in its neutral and ionized states, posits that internal modes of the 
olfactant are excited impulsively during hole transfer from a donor to acceptor site on the receptor, specifically 
those modes that are resonant with the tunneling gap. By projecting the impulsive force onto the internal 
modes, we can determine which modes are excited at a given value of the donor-acceptor tunneling gap. Only 
those modes resonant with the tunneling gap and are impulsively excited will give a significant contribution to 
the inelastic transfer rate. Using acetophenone as a test case, our model and experiments on D. melanogaster 
suggest that isotopomers of a given olfactant give rise to different odorant qualities. These results support 
the notion that inelastic scattering effects play a role in discriminating between isotopomers, but that this is 
not a general spectroscopic effect. 



I. INTRODUCTION 

The general model for detection is that the response 
is triggered by the transfer of an electron from a donor 
(D) to an acceptor (A) within the receptor site by the 
presence of an olfactant molecule that provides a bridge 
between the two. In the absence of the odorant, the dis- 
tance between D and A is too great and electron transfer 
is inefficient. Placing an olfactant (B) between the two 
allows the electron transfer to occur either as a single 
coherent scattering event from D to A, or as a sequence 
of two incoherent hops, first from D to B then from B 
to A. This simple model is consistent with the standard 
"swipe-card" model for odor detection since implicit in 
this is that B must fit into some sort of pocket and be 
in a correct alignment between D and A in order for the 
charge transfer process to occur. 

An interesting extension to this paradigm is that 
molecular shape may not be the sole deciding factor in 
scent recognition. Indeed, it has been observed recently 
that fruit flies {D. melanogaster) can naively discriminate 
between several isotopomers that have deuteriums sub- 
stituted for hydrogen^. These flies can also be trained 
to associate a specific odorant with an electric shock 
resulting in the specific avoidance of the conditioned 
odorant^. Flies trained in this manner can discrimi- 
nate deuterated from hydrogenated isotopomers, indicat- 
ing an ability to perceive differences in these odorant s. 
Mutant flies that lack the orco ofactory co-receptor gene 
are broadly anosmic and fail to discriminate between iso- 
topomers, indicating the behavioral discrimination of iso- 
topomers relies on olfactory perception. Lastly, using 
the negatively- reinforced learning paradigm, Drosophila 
has been shown capable of generalizing a trained odor- 
ant to another that shares a similar vibrational spec- 
trum. For example, when deuterated d5-benzaldehyde 
is paired with electric shocks, the fiies will subsequently 



avoid fully-deuterated dl7-octanol (vs. undeuterated h- 
octanol), suggesting that the C-D vibrational mode is a 
salient feature of odorant perception in this species. 

A controversial and speculative explanation of these 
observations by Turin is that there is an spectroscopic 
component to olfaction. ^'^ Turin's model supposes that 
the bridge facilitates inelastic electron scattering and as 
such the olfactory response can be predicted by com- 
paring infra-red (IR) spectra of various olfactants. Iso- 
topic substitution of H for D shifts the CH stretching 
frequency from the 2850 - 3100 cm~^ range into the 2300 
cm~^ range. Few molecules absorb in this IR region and 
there is little or no biological need or evolutionary pres- 
sure that we know for detecting deuterated compounds. 
Needless to say, this theory has been met with consider- 
able skepticism since it runs contrary to more commonly 
held model in which the geometric shape and chemical 
nature of the olfactant are the primary components of ol- 
factory reception.^ Moreover, many of the claims in Ref. 
3 and Ref. 4 were show to be inconsistent with psy- 
chophysical tests performed on human subjects by Keller 
and Vosshall.^ For example, in this study human subjects 
could not discriminate by smell between deuterated and 
non-deuterated forms of acetophenone. 

While the olfactory receptors and the co-receptor have 
been identified and is believed to be a multi-helix trans- 
membrane protein with the binding sites on the periplas- 
mic domain and a G protein-binding site site on the cy- 
toplasmic domain, information concerning the molecular 
structure at the olfactant binding site remain illusive, but 
appear to involve odorants binding to Cu(II) and Zn(II) 
ions bound to the protein loop.*" We can assume that 
at the heart of the process, a charge is transferred from 
a donor to an acceptor within the receptor site itself. 
In the absence of an odorant, the distance between the 
donor and acceptor is too great to allow efficient elec- 
tron transfer. When the odorant is present, it acts like a 
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bridge between the donor and acceptor allowing for more 
efficient electron tunneling between the two.^ Aromatic 
compounds are typically better electron donors {i.e. hole 
receptors) due to their electron-rich tt -system. With this 
in mind, we shall assume that the bridge/olfactant acts 
as an intermediate for hole transfer between a donor site 
and an acceptor site. Energy conservation requires that 
the energy transferred to vibrational motion is equal the 
tunneling gap. 



II. THEORETICAL MODEL 

For a donor-bridge-acceptor system, the tunneling ma- 
trix element for electron transfer rate can obtained from 
second-order perturbation theory^'^ 

VdaiEda) = {A\VAB{Eda " Hb)-^Vbd\D) (1) 

where \A) and \D) denote the final and initial states and 
Hb is the electronic Hamiltonian for the bridge molecule. 
Here, an electron can tunnel from an occupied molecular 
orbital on the donor, to an unoccupied molecular orbital 
on the bridge, and finally be transferred to an unoccupied 
molecular orbital on the acceptor, E^a is the tunneling 
gap and Vij denotes the couplings between the donor to 
bridge and bridge to acceptor. One can also consider hole 
transfer along similar lines. 

We define the couplings between the donor to bridge 
and bridge to acceptor as Vdb = Jdb\d){b\ and V^a = 
Jba\b){a\^ respectively. From the Franck-Condon prin- 
ciple, these can be factored from the nuclear/vibrational 
terms. Further, let us assume that the bridge has two 
electronic states, denoted as |1) corresponding to the neu- 
tral bridge, and |2) corresponding to the oxidized bridge. 
In addition, the bridge has a set of internal molecular vi- 
brational states local to the molecule with normal mode 



frequencies, {ujq}. It is important to point out that the 
vibrational states of the neutral are not the same as those 
of the ion due to shift of the energy minima and (Duschin- 
sky) rotation of the normal modes. 

We can thus evaluate the transition matrix as 



Vda = JdbJba ^ 



(n;|n2)(n2|0i) 



n2,n'^ 



Eda - {€2 ^huJ2 -{112 ^1/2))'' 



(2) 



where the sum is over all intermediate, n2, and final, 
n'^, vibrational states (including overtones and combi- 
nations) of the oxidized bridge molecule and we as- 
sume that the initial vibrational populations are in their 
ground states. The terms in the numerator are the 
overlap integrals between the vibrational states of the 
neutral and oxidized bridge. In general, these overlaps 
can be computed using the approach outlined in Sharp 
and Rosenstock^^'^^; however, this formally this involves 
summing over all possible combinations of initial, final, 
and intermediate vibrational quantum numbers including 
overtone modes. 

To simplify this, let us consider the bridge ioniza- 
tion/deionization process to be purely impulsive so far 
as the nuclear dynamics are concerned. Upon ioniza- 
tion, the nuclei in B experience a sudden change in their 
electronic environment. Ordinarily, the vibrational wave 
function would begin to evolve on the Born-Oppenheimer 
potential, V^~^\ corresponding to this new electronic con- 
figuration. However, before nuclear motion can occur, 
the charge hops to the acceptor returning B to its neutral 
electronic state. Over the course of the charge-transfer 
event, the nuclei in the bridge molecule receive an impul- 
sive "kick" due to the sudden changes in electronic state. 
This excitation can be imparted to the nuclear motion 
on the ground-state (B) potential by an instantaneous 
shift in the momentum of the ground-state vibrational 
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FIG. 2. Olfactory Discrimination of Acetophenone isomers, a) Drosophila melanogaster (it;^^^^) naively avoid ACP in a T-maze. 
The behavioral responses to both d5-ACP and d8-ACP are significantly stronger than to both the d3-ACP and h-ACP. (n = 22 
groups each) (b) Drosophila can discriminate between d5-ACP and d3-ACP after training with electric shock. During the 
training session, populations of flies are exposed to one of the isotopomers paired with an electric shock. When subsequently 
tested in the T-maze, the flies significantly prefer the unpaired odorant over the paired odorant {p < 0.0001; n = 44 groups 
each), (c) Drosophila can discriminate between h-ACP and d3-ACP after training with electric shock. 



wavefunction along the direction of VV^~^\ 

y.. = %^5:(n;|e^nOi). (3) 

It is important to recognize the p term in this expression 
corresponds to the momentum imparted to the nuclear 
motion due to the transition and that energy conserva- 
tion requires E^a = /2m be the additional kinetic en- 
ergy given to the internal nuclear motion of the olfactant 
bridge molecule. We can factor this into contributions 
from each normal mode and require that the momentum 
imparted to a given vibrational mode is proportion to the 
energy gradient along that direction. 

The immediate advantage of this model is that we only 
need to consider the nuclear motion of the neutral species 
and in practice we can evaluate the matrix element ex- 
actly. However, it is useful to examine just the leading 
order terms. If the momentum kick is small, then we can 
expand the translation operator for each normal mode as 

(n'iJe*P-->i,) = 6on + ip- (n'ljf |0i,). (4) 

This first term is simply the elastic component, which 
is expected to dominate. The second term contains the 
matrix element of the position operator. This implies 
that to a first approximation the inelastic scattering of 
an electron via the bridge species excites the vibrational 
modes of the bridge that are also infra-red (IR) active. 
This is a central component of Turin's spectroscopic the- 
ory of olfaction. ^'^'^ However, only those modes that are 
directed along VF^^^ will be excited by the impulsive 
scattering process. 

Equation (4) gives us a direct way to rapidly screen 
whether or not a given odorant is expected to exhibit 
an isotope effect using quantum chemical means. Start- 
ing from the equilibrium position of the neutral, one first 



determines the vibrational normal modes and frequen- 
cies of a given olfactant. This gives the IR response of 
the molecule. We then determine the energy gradient 
of the ionized species at this geometry and project this 
onto the normal modes. Since E^a remains an unknown 
in our model, we simply require that this energy be dis- 
tributed amongst the normal coordinates in proportion 
to their projection onto the energy gradient. Thus, even 
if a mode has a strong IR response and satisfies energy 
conservation, unless that mode is directly excited by the 
impulsive scattering process, then it will not contribute 
to the overall transition rate. Our approach is similar to 
that of Brookes et al}'^] however, our approach explicitly 
treats the olfactant as a bridging species between the 
donor and acceptor and derives all relevant input from 
quantum chemical treatments of the bridge. This allows 
us to make specific predictions regarding the potential 
olfactory response of a given scent molecule. 

As a test case, we consider acetophenone (ACP) and 
its deuterated isotopomers. In d3-ACP, the hydrogens 
on the methyl group are replaced, in d5-ACP the hy- 
drogens on the phenyl ring are replaced, and in d8-ACP 
all hydrogens are replaced with deuteriums. In all cases, 
we first perform geometry optimizations of the neutral 
species followed by calculations of vibrational frequencies 
using the B3LYP density functional with 6-31 G(d) basis 
set in vacuum using the NWChem quantum chemistry 
package. The same functional and basis set were then 
used to calculate the nuclear gradients and frequencies 
of the radical cation at the equilibrium geometry of the 
neutral. This approach is numerically reliable and pro- 
duces accurate potentials and gradients for this system 
without undue effort. Moreover, this general procedure 
can be applied to a broad set of olfactants. 

In Figure 1(b) we show the predicted IR response and 
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the projection of the energy gradient of the oxidized 
species onto the normal mode frequencies of the neu- 
tral. Briefly, the peaks around 3100 cm~^ correspond 
to the C-H stretching modes of the molecule. In the fully 
deuterated d8-ACP, this peak is shifted to 2300 cm~^, 
consistent with the typical isotope shift of a C-D stretch- 
ing mode. Furthermore, the IR spectra for all but the 
d8-ACP is void of peaks between 1800 - 3000 cm-^ Only 
the fully deuterated d8-ACP has any appreciable spectral 
density in this region. While there is little to note in com- 
paring the intensities for the modes below 1800 cm~-^, the 
contribution to the gradient coming from the C-H and C- 
D stretching modes is remarkably strong. Most notable, 
the peak at 2300 cm~^ is more or less the sum of contri- 
butions from the ring and methyl C-D stretching modes. 

Based upon the relative intensities of both the IR and 
energy gradients around 2300 cm~^ in Figure 1(b), the 
electron transfer model suggests that D. melanogaster 
should readily discriminate h-ACP and d8-ACP. In com- 
paring the predicted IR spectra, there is almost no IR 
oscillator strength in the C-D stretching modes for d3- 
and d5-ACP while d8-ACP exhibits a fairly strong sig- 
nal. However, in comparing the gradients, the d5-ACP 
gradient has considerably smaller projection onto the C- 
D stretches modes and the d3-ACP projection is almost 
vanishingly small. Based upon this, one predicts that 
there should be an ability to differentiate between h-ACP 
and dS-ACP. Moreover, the flies may be able to discrim- 
inate between d3-ACP and d5-ACP even though both 
have very weak IR signals at 2300cm~-'^. 



III. EXPERIMENTAL TESTS 

Previously, d5-ACP and d8-ACP odorants were found 
to be signiflcantly more aversive to Drosophila than h- 
ACP^. We verified this result, examining the naive avoid- 
ance of Drosophila w^^^^ to 0.3% h, d3, d5, and d8 ACP 
using an olfactory T-Maze (Figure 2a). The avoidance 
response to h-ACP and d3-ACP were indistinguishable 
(F3,80 = 11.724, Bonferonni-Dunn post hoc p = 0.276); 
the responses to d5 and d8 were also not significantly 
different from each other (p = 0.842). However, the 
avoidance response to d5-ACP was significantly greater 
than the response to d3-ACP {p = 0.0016), suggest- 
ing deuteriums on the benzene ring provided a different 
odorant quality than the methyl deuteriums, leading to 
an increased saliency for the d5-ACP and d8-ACP. To 
test this hypothesis, we used negatively-reinforced olfac- 
tory learning^"^. Surprisingly, the w^^^^ flies are capa- 
ble of discriminating between d5-ACP and d3-ACP after 
training, indicating a signiflcant difference in perceptual 
quality between these odorants (Figure 2b; t = 24.714, 
p < 0.0001). The w^^^^ flies are also capable of dis- 
criminating between h-ACP and d3-ACP after training 
(Figure 3c; t = 10.596; p < 0.001). The conditioned 
avoidance in this experiment was reduced in comparison 
to the conditioned avoidance of d5-ACP vs. d3-ACP, in- 



dicating that d3-ACP is more similar in odorant quality 
to h-ACP than to d5-ACP. The robust discrimination 
between these very closely related odorants is consistent 
with our model that predicts differences in the excitation 
of vibrational modes of the methyl and deuterium C-D 
bonds by the impulsive scatting. 



IV. SUMMARY 

It is clearly evident that D. melanogaster can be 
trained to distinguish deuterated isotopomers of vari- 
ous olfactants. Whether or not this is the result of a 
spectroscopic detection mechanism is an another ques- 
tion. While deuteration does shift the C-H stretch- 
ing mode into the otherwise empty 2300 cm~^ region, 
our model suggests that IR response alone does not de- 
termine whether or not the flies can discriminate iso- 
topomers. Moreover, this effect appears to be limited to 
a single spectral region and does not involve the entire IR 
spectrum as suggested by Turin. ^'^ A crucial test of our 
model would be to identify an olfactant that can not be 
discriminated from its deuterated isotopomers. For ex- 
ample, in removing an electron from ethylene, the C=C 
bond-length in C2H4 will be longer than in C2H4 (due 
to reduction in bond-order) but the C-H bond-lengths 
will be more or less unchanged. This suggests that the 
energy gradient in our model will be directed primarily 
along the C=C bond and have very little projection onto 
the C-H stretching modes. However, the IR spectrum of 
ethylene does show an isotope shift of the C-H stretch- 
ing modes into the 2300cm~^ region. Hence, we predict 
that D. melanogaster should not be able to discriminate 
between ethylene and its isotopomers. 

Lacking detailed molecular level knowledge of the re- 
ceptor sites, we speculate the tunneling gap between 
donor and acceptor sites is around 2300 cm~^ and that 
the narrow spectral region about E^a is giving rise to the 
observed ability for Drosophela to discriminate between 
deuterated isotopomers by scent. It is possible that the 
tunneling gap may be tuned to the nitrile (-C=N) stretch- 
ing mode which does lie in this spectral region and may 
serve a role in detecting a variety of naturally occurring 
nitrile containing odorants. 
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